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ABSTRACT 

The accuracy of ground-based astronomical photometry is limited by two factors: pho- 
ton statistics and stellar scintillation arising when star light passes through Earth's 
atmosphere. This paper examines the theoretical role of the outer scale Lq of the op- 
tical turbulence (OT) which suppresses the low-frequency component of scintillation. 
It is shown that for typical values of Lq ~ 25 — 50 m, this effect becomes noticeable 
for a telescopes of diameter around 4 m. On extremely large, 30 — 40 m, telescopes 
with exposures longer than a few seconds, the inclusion of the outer scale in the 
calculation reduces the scintillation power by more than a factor of 10 relative to con- 
ventional estimates. The details of this phenomenon are discussed for various models 
of non-Kolmogorov turbulence. Also, a quantitative description of the influence of the 
telescope central obscuration on the measured scintillation noise is introduced and 
combined with the effect of the outer scale. Evaluation of the scintillation noise on the 
future TMT and E-ELT telescopes, predicts an amplitude of approximately 10 /imag 
for a 60 s exposures. 
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1 INTRODUCTION 

Stellar scintillation is the random fluctuation of the radia- 
tion flux entering the aperture of the telescope caused by am- 
plitude distortions of light wave passing through the turbu - 
lent terrestrial atmosphere {Tatarskii 1967: [Roddier|[T98ll ). 
The photometric error due to scintillation, th e scintillation 
noise, has been rep eatedly studied (see, e.g., lYound [19691 : 
iDravins et al.| [l998) , as it often determines the fundamental 
limit of the accuracy of ground-based photometry. 

The basic dependencies required to calculate scintil- 
lation noise have been known by th e wide astronomical 
community for qui te a while now l|Heaslev et al. 1996; 
iGilliland et al. Ill 9931 ). However, these relationships were ob- 
tained for ideal cases and cannot always be used for accurate 
prediction of the scintillation noise. This became especially 
noticeable when accurate measurements of the intensity of 
the optical turbulence (OT) in the atmosphere at altitudes 
responsible for the occurrence of scintillation on la rge tele- 
scopes, began to be accessible l|Kenvon et al.ll2006t ). 

The content of this paper is a theoretical study of two 
factors affecting the power of stellar scintillation on large 
and extremely large telescopes: the influence of the OT outer 
scale and the telescope central obscuration (CO). 

Section [5] recalls the basic description of the phe- 
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nomenon and its relationship with several parameters of the 
OT in the atmosphere. In the following section, we assess 
the influence of the outer scale in the case of different and 
simplified non-Kolmogorov models. The effect of the CO is 
considered in Section 3] first for the Kolmogorov model, and 
then for a general case. In the last Section, the conclusions 
are formulated and a prediction of the scintillation noise for 
future extremely large telescopes is given. 



2 THEORY 

2.1 Scintillation noise in photometry 

Measurements of the brightness of astronomical objects are 
always burdened by noises of several different nature. De- 
pending on the origin, these noises are included in the mea- 
sured signal in different ways. Noise caused by the stellar 
scintillation, enters by multiplicative way, i.e. for it, the sig- 
nal to noise ratio is independent of brightness. Stellar scin- 
tillation is usually characterised by index s 2 , representing 
the variance of the relative fluctuations of the intensity / of 
the radiation, passing through the receiving aperture: 



<(/-(/» 2 >/</r 



(i) 



Although in this definition, the averaging over all possible 
states (over ensemble) is performed, in practice, using the 
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ergodic property for the turbulent phenomena (jTatarskiil 
Il967f ). it is replaced by averaging over time. 

For astronomical applications, the telescope entrance 
pupil D is considered as a receiving aperture. The scintil- 
lation noise can be expressed in magnitudes a 2 ^. It then 
becomes the additive noise, associated with the scintillation 
index: — 1.179 s 2 . 

In nighttime conditions, suitable for astronomical obser- 
vations, the OT is usually such that the phase and amplitude 
distortions at the pupil plane are well described und er the 
approximation of weak perturbations (Tatarskii 1967). It is 
believed that for large apertures (geometric optics regime), 
this is even more so. In this approximation, the turbulence 
of each layer is independent of the previous layers, and the 
scintillation index observed on the surface, is linearly re- 
lated to the distribution of the structural coefficient C 2 (z) 
on t he line of sigh t throughout the whole atmosphere (see, 
e.g.. lRoddierlll981f ) by the following integral: 

s 2 = / C 2 (z)Q(z)dz, (2) 

J A 

where z is the distance to the turbulent layer, which in the 
case of measurements at zenith, coincides with the altitude h 
above the observatory. The weighting function Q(z) relates 
the output from the layer to its corresponding effect at the 
surface. 



2.2 Basic relations 

The calculation of t he weighting function Q (z) has been de- 
scribed in detail in (|Tokovininl 120021 . 120031 ). It involves the 
integration of the 2D spatial spectrum of the amplitude per- 
turbation over all frequencies. In the case of axial symmetry 
of integrand functions, it is easier to perform the integra- 
tion in polar coordinates, where, after averaging over polar 
angle, the functions depend only on the modulus / of the 
vector of spatial frequencies. For isotropic and locally homo- 
geneous OT with a spatial spectrum of the wavefront phase 
fluctuations $(/), normalized to the value of C 2 (z): 

roc 

Q(z) = 9.61 $(f)S(z,f)A(f)fdf, (3) 
Jo 

where the Fresnel filter S(z, /) describes the evolution of 
the amplitude perturbations in the propagation of the light 
wave. In the case of monochromatic radiation with wave- 
length A, the filter S(z,f) — sin 2 (irXzf 2 )/X 2 , and its intrin- 
sic spatial scale is defined by the Fresnel radius tf = (Xz) 1 ^ 2 . 

Aperture filter A(f) takes into account the averaging 
of wave amplitude by a receiver. For a circular entrance 
aperture D it is the Airy function (2Ji(ivDf)/(iTDf)) 2 if 
there are no other factors for averaging. However, in general, 
stellar scintillation are registered as temporary fluctuations 
in light intensity averaged over a measurement (exposure). 
Dependence of the scinti llation on th e exposure was inves- 
tigated repeatedly (e.g., lYound Il967l ; iDravins et al.l I1998T ) 
with the help of time averaging the signal. 

On the basis of the "frozen turbulence" hypothesis 
|Tavlorlll93sl ). the temporal averaging can be replaced with 
a spatial filt e ring of the scint i llation spectrum (see, e.g., 
iMartinl 1 19871 : iTokovininl 120021 : iKornilovl |201lD . extending 
the concept of aperture filtering. The scintillation power in 



the case of a finite exposure is still described by the ex- 
pression ([3]) if instead of A(f) we substitute t he product 
A(f) A s (f), where A s (f) is the wind shear filter (|Tokovininl 
l2002l : lKorniiovll201ll '). The real "freeze" of the OT is not re- 
quired as the invariability of the spatial spectrum is enough. 

In the cases where measurements are performed with an 
exposure r so short that the wind w in the atmosphere, does 
not shift the OT by a significant distance tw <C D, the addi- 
tional averaging can be neglected (filter A s (f) is much wider 
than the filter A(f) and multiplication does not change the 
integrand). Hereafter, this situation is called short (zero) 
exposure regime (SE). 

In real astronomical observations, the opposite situa- 
tion is much more common: during the exposure, the wind 
shifts the OT by distances exceeding the aperture of the 
telescope tw 3> D, and temporal averaging becomes the 
dominant effect (filter A 3 (f) is narrower than A(f) and it 
defines the integrand). This case is denoted below as long 
exposure regime (LE). 

2.3 Large apertures approximation 

The expression Q is valid for arbitrary parameters, but in a 
general form it can be integrated only numerically. To study 
the dependence of Q(z) on distance z and other parameters, 
it is necessary to perform some simplification. 

To describe the scintillation on a typical telescope, the 
approximation of large aperture J) > rF can be used. This 
approximation is well satisfied for telescopes with diameters 
of 1 m or more, since for any reasonable atmospheric condi- 
tion, tf < 0.1 m in the optical and near-IR spectral range. 

In this case, one can use a well-known feature of the 
spectral filters included in the integral (|3]): the filters A(f) 
and S(f) overlap a little. After replacing sin(7r Az/ 2 ) on 
hvXzf 2 ), the integrand is still largely unchanged l|Roddierl 
1981). As a consequence, the dependence on A disappears 
not only for monochromatic case, but in the case of a wide 
spectral band of d etector, so the m ore accurate description 
of this filter from ||Tokovininl l20031 is not necessary. 

Passing to the dimensionless frequency q = fD, we can 
write for circular aperture and the SE regime: 

r oo 

Q s (z) = 38A4D- 7/3 z 2 4>(q) q 3 ( J 1 (7rg)) 2 dq, (4) 
Jo 

In the LE regime, the approximation technique is par- 
ticularly suitable because the wind smoothing additionally 
suppresses the high-frequencies. We use the asymptotic be- 
havior for the c orresponding f ilter of the wind shear A s — 
D/nrwq from (|Kornilovll201ll ). Multiplying the integrand 
by it, we obtain: 

r oo 

Ql(z) = 12.24 D~ 4/3 z 2 (rwy 1 / $(g) q 2 ( Jj (7vq)) 2 dq, 

Jo 

(5) 

It should be noted again that the conventional depen- 
dencies on telescope diameter D and propagation distance 
z are the result of the the Fresnel filter approximation un- 
der the condition D ^> r-p and are not associated with a 
particular form of the OT spectrum. 

For the Kolmogorov model of normalized spatial spec- 
trum of phase perturbations &(q) = g _11,/3 , the inte- 
grals Is and II in the right-hand side of these expres- 
sions, are easily calculated and are equal to 0.4508 and 
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0.8699, respectively (see equations IA1I and IA2|) . As a re- 
sult, the previous expressions reduce to the known depen- 
dencies Qs(z) = 17.34 D~ 7/3 z 2 for SE regime and Q L (z) = 
10.66 D~ 4 / 3 z 2 (tw)~ 1 for LE regime. After integration over 
the whole atmosphere, they le ad to the fol l owing formulae 
of observed scintillation noise l|Youndll967l : iGilliland et alJ 
ll993l ; lKenvon et alj|2006h : 

4 = 17.34 D~ V ' A [ d(z)z 2 dz, (6) 

J A 

and: 

2 i n cc 71-4/3 -1 f Cn( z ) Z 2 j /»\ 

s L = 10.66 D ' r / — dz. (7) 

Ja h>[z) 

The integrals in these formulae (the atmospheric moments) 
are determined by the particular state of the atmosphere 
and are the figures of merit wh en monitoring and /or forecast 
of co nditions for photometry l)Kenvon et al.ll2006l : iKornilovl 
The matter of this paper is to consider the effects 
which modify the coefficients before the atmospheric mo- 
ments. 



3 THE IMPACT OF THE OUTER SCALE OF 
TURBULENCE 

The immediate perception of the scintillation as a result of 
observations with a naked eye and numerous experiments, 
performed with small telescopes, have created the stable il- 
lusion that the scintillation is a sufficiently high frequency 
process. The problem of contribution of the high- and low- 
frequencies in th e scintillation has been discussed for a long 
period (see, e.g.. lYoungl fl967'). E.g., in one of the first pa- 
pers devoted to stellar scintillations l|Reiger|[l963l '). the au- 
thor claimed that turbulence outer scale has no effect on the 
scintillation intensity, in contrast to the inner scale. The fact 
that the situation changes radically when the measurements 
are made on large telescopes, usually goes unnoticed. 

Nevertheless, we decided to re-evaluate their relative 
contribution. Assuming the Kolmogorov spectrum, the frac- 
tion of the scintillation power in the frequency interval be- 
tween and some dimensionless frequency q was calculated. 
The calculation shows that most of the scintillation power 
(~ 70%) passes through the first (main) peak of the filter 
A(q). In the LE regime, the fraction of the power for this 
peak is even greater (~ 95%). The main peak of the aperture 
filter is located at the q < 1.22, which for large telescopes 
corresponds to meters scale. It is logical to expect that the 
distinction of the real OT from the usual Kolmogorov model 
at such scales, may significantly affect the measured power. 

3.1 Non-Kolmogorov models with outer scale 

Skipping the discussion of the physical meaning of the outer 
scale of turbulence Lq l|Tatars kii 1967), which plays an im- 
portant role in the generation of turbulence per se, we will 
consider it simply as an additional parameter in the mathe- 
matical description of the spectrum of perturbations in the 
refractive index. The main purpose of this parameter is to 
limit the infinite spectral density at frequency / = inherent 
in the Kolmogorov spectrum, or to overcome the divergence 
in the relationships derived from this spectrum. 
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Figure 1. Spectral power density of the phase perturbations as 
the function of dimensionless frequency q for different models of 
non-Kolmogorov turbulence for different values of the dimension- 
less outer scale v. Solid lines depict the VK model, short dashed 
the GT model, long dashed the EM model, dash-dot lines indi- 
cate the ME model. See abbreviations in the text. (Online version 
colored) 

In astronomical applications, the outer scale of the 
OT was studied for a l ong time in connection with long- 
baseline interferom etry (|Davis et al.lll995l ; lAvila et al ] |l997l ; 
iMaire et al.l 120061) . The development of adaptive optics on 
large telescopes also requires the correct description of the 
spect rum of phase perturbations in the low-frequency re- 
gion dConan et al.ll2003l ; iTokovinin. Sarazin fc SmettefeoOTj ; 
iMartinez et al.ll201Ch . 

The most generally accepted OT model w hich includes 
the outer scale is the V on Karman (VK) model |von Karmanl 
ll948l;lTatarskiilll967l ). In this model, the normalized spec- 
trum of phase perturbation can be written as: 

4>(/,£o) = (/ 2 + Lo 2 r 11/6 - (8) 

The main feature of the VK model is a saturation of the 
power density in the range / < Lq 1 at the level of Lj 1 ^ 3 . 

An alternative model without saturation but with a 
change in exponent of po wer law, was proposed by Gree n- 
wood and Tarazano (GT) l|Greenwood k, Tarazanolll974h : 

«f(/,io) = (/ 2 + /L - 1 )- 11 / 6 . (9) 

It is easy to see that at low-frequencies, the spectrum has 
an asymptote ~ Lj 1//6 / _11//6 . The intermediate range is very 
wide, so in this model, the divergence from the Kolmogorov 
power spectrum become apparent early enough. 

The expo nential model (EM) beh aves similarly, it was 
considered in (|Lukin fc Pokasov| [l981) as follows: 

4>(/,L ) = r 11/3 [l-cxp(-/ 2 L ( 2 ) )]- (10) 

This spectrum is asymptotic to ~ Lq/ _5//3 at low frequen- 
cies, ie, has a slightly smaller slope than the GT models. 
The intermediate range of the spectrum is quite narrow. 

We modified the EM model so, that the spectral density 
at the origin is finite and coincides with ij 1//3 , predicted by 
the VK model. This exponential model (ME) is described 
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by the dependence 



Hf,L ) = r 11/3 [i 



exp(-/ L 



11/6 ,-11/6^2 



(11) 



The only difference between this spectrum and VK model is 
a faster transition from one asymptotic branch to another. 

The behavior of the spectral densities for these four 
models are shown in Fig. [T] as a function of dimensionless 
frequency q = fD for different dimensionless outer scale 
frequency v — D/Lq. The vertical lines indicates the fre- 
quency domain where the main peak of the aperture filter is 
located. Note that for any model, the spectral density can 
be described as a product of Kolmogorov spectrum by some 
spectral filter suppressing low frequencies. 

The question of which model best describes the 
OT in the atmosphere is still open, despite many ex- 
perimental attempts to solve it (e.g., iMaire et al.l 20081 ; 
IWheelon. Short fc Townesll2007r ). For this reason, in further 
calculations, we use all of the above model with sufficient 
variety of behavior. 



3.2 The effect of the outer scale on scintillation 
power 

The impact of the outer scale can be achieved by substi- 
tuting the appropriate model in the formulae (f4]) and ([5]) 
expressed as a function of dimensionless frequency q — fD 
and outer scale v = D/Lq. It is obvious that after such 
replacement as well as in the case of the Kolmogorov spec- 
trum, dependence on the parameters z and D can be taken 
out of the integral completely. In the case of the SE regime 
we obtain: 



i(z,u) = 38.44 D 



- 7/ v 



g 3 (Ji(7rg)) 2 dg. (12) 



A similar expression can be written for Ql(z, v) in the case 
of LE regime. The value of the integral in the formula de- 
pends explicitly only on the parameter v. However, if the 
outer scale Lo varies with altitude, then there is an implicit 
dependence of the integral on the altitude and the diameter 
of telescope via the parameter v. 

The approximation for large aperture ((4l [5j can be used 
if not only the D ^> tf but also Lo 3> rw. Otherwise, the 
integrals diverge and a formula with the exact expression for 
the Fresnel filter should be used. In that case the integrand 
maximum is located near ~ 1 /r? and the situation becomes 
similar to the scintillation in a small aperture. 

Multiplying and dividing the right side of (|12l) by the 
value of the integral Ts we get that Q(z,v) = Q(z)G(u), 
where the function G(y) describes the impact of the outer 
scale, and in SE regime is equal to 



Qs{v) =Z S 



$(g,/.)g 3 (Ji(^)) 2 dg. 



(13) 



The smaller this term is, the greater is the effect. The possi- 
ble factorization of the weighting functions leads to the fact 
that the scintillation index s 2 for non-Kolmogorov OT con- 
nects with the index s 2 , defined by the formulas © or ([7]) 
by the simple relationship: s 2 = s 2 Q(v). 

In principle, for exponential models the integral ()13[) 
can be calculated analytically, but the resulting expression 
is so cumbersome that it is meaningless. For VK and GT 
models we have failed to get an analytic solution, and the 




Figure 2. Dependences of the integrals Q (y) (thin solid and 
dashed lines) and Q+(v) (dot-dashed line) and the sum of Q(y) 
(thick lines) on the dimensionless frequency v for the piecewise 
model. Solid lines depict the saturated model, the dased the non- 
saturated model. On the left: SE regime, on the right: LE regime. 



usage of expansions of the type jMaire et al.ll2007l ) leads to 
the appearance of improper integrals. 

We therefore investigated an overall behavior of Q(y) us- 
ing the piecewise power-law dependence with a break point 
at a frequency v. The necessary integrals are calculated an- 
alytically within two segments of {0, v} and {v, oo}. We will 
refer to the integral with the spectrum of g -11 ^ 3 in the range 
{is, oo} as Q + {v) after its normalization with I. The inte- 
gral with the spectrum i/ _11//3 (the model with saturation) 
or v~ 2 q~ 5 ^ z (the model without saturation) in the range 
{0, v} will be referred to as Q~ (y). Those constants provide 
continuity to the spectrum at q = v. We will specify the SE 
or LE regimes, using the appropriate subscript. The integra- 
tion results of these functions are given in Appendix I A2I 

The functions G~(v) and G + (v) are shown in Fig. [2]for 
both regimes. It is evident that the decrease in scintillation 
power caused by the outer scale, arises from a sharp fall in 
the function Q + (u), which is not compensated by an increase 
in Q~(y). The figure also shows their sum, i.e. dependence 
of Q(y). At the point v ~ 0.5 the power falls at w 20%. 
As expected, the effect in the case of saturated spectrum is 
slightly more than for the non-saturated one. 

The asymptotes of the Qs (v) for small v follows from the 
evaluated integrals in Appendix IA2I and properly describes 
the functions while v < 0.4: 



g a (v) = 1 - 1.433 u v ' A + 1.428 z/ J/d + 



13/3 



9s{y) 



1.083 v lz +0.984 v 13/3 



+ ■ 



<h 



-2 -5/3 

v q ' 



(14) 



For large u, these functions come to the asymptotes 
0.412 v~ 2/3 and 0.508 ^" 2/3 , respectively. 

In the LE regime, the effect of the outer scale is much 
stronger: at u w 0.5 the power is reduced by almost half and 
at v ss 1 about 5 times. The approximation of the functions 
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Figure 3. Dependence of the Q(y) on the dimensionless outer 
frequency v for different models of non-Kolmogorov turbulence 
for SE regime. Solid line depcts the VK model, short dashed for 
the GT model, long dashed for the EM model, dash-dot lines 
indicates the ME model. 



Ghiv) looks as: 

g h (v) = 1 - 1.560 u 4/s + 1.099 v 10/3 + ..., $ = v - lirA 

Q^{v) = 1 - 1.276 u 4/3 + 0.787 u 10/z + ..., $ = v~ V 5/3 , 

(15) 

and approximates them into region v < 0.4 with an accuracy 
better than 0.02. In the v > 1, Ql{v) tend to 0.128 u~ 5/3 
and 0.375 u~ 5 ^ 3 , respectively. 

The difference between the behavior of the functions 
evaluated for the spectra with and without saturation is not 
fundamental and can hardly be detected in actual measure- 
ments. As we are to see from calculations using the models 
described in section [3"TTl a much greater effect occurs because 
of the length of intermediate region in these models. 

Fig.[3]and Fig.[4]show the functions G{y), calculated by 
numerical integration for different non-Kolmogorov models. 
It is seen that the initial parts (y < 0.2) of the curves differ 
greatly in both SE and LE regimes. E.g., at the point v = 
0.1, the effect is almost zero for the EM model but for the GT 
model, which has the longest transition region, the function 
Qiy) « 0.75. In the LE regine, the Q(y) is 0.92 and 0.60 
respectively for these two models. Between these extremes, 
the curves pass for VK and ME models. 

The absolute difference between the models decreases 
with increasing v. For models with saturation, the curves 
pass below the others as seen in Fig. 2] A comparison with 
Fig. [2] shows that the piecewise power-law model approx- 
imates a general shape and, having zero transition range, 
can be considered as a lower estimate of the effect (or up- 
per bound for the Q(y) function) for an appropriate type of 
models. 



4 THE IMPACT OF THE CENTRAL 
OBSCURATION 

Most astronomical telescopes have a central obstruction 
(CO) caused by the secondary mirror, therefore its im- 
pact on power and spectrum of the scintillation was repeat- 
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Figure 4. Dependence of the Q{v) on the dimensionless outer 
frequency v for different models of non-Kolmogorov turbulence 
for LE regime. See designations in Fig. [3] 



edly considered (|Youndll967 l: lDravins et alJ l998). However, 
these studies have been focused more on qualitative assess- 
ment. Calculations of the scintillation spectra served mainly 
to showcase the presence of the effect. There are much more 
studies of the effect of various occul ting structures on diffrac - 
tion image in telescope focal plane l|Harvev fc Fta clas 1995), 
which mathematically coincides with the aperture filter nec- 
essary for us. 

In this section, we once again examine this effect in 
the approximation of large telescopes, and obtain the ex- 
pression suitable for accurate description in the case of the 
Kolmogorov spectrum. 

4.1 Central obscuration for Kolmogorov spectrum 

In the presence of the CO, the aperture filter A(f) 
is described by the well-known expression (see, e.g., 
lHarvev fc Ftaclasl Il995l ) depending on additional parame- 
ter e, which is a ratio of inner diameter of entrance pupil 
to its outer diameter. In the formula Q and follow-up, the 
aperture filter is separated into factors, and the factor 4/7T 2 
is taken out of the integrals. Therefore, we write only the 
part a(q,e) of the filter, which should replace the (Ji(7rq)) 2 
in the integrands. 

a(g.e) = T(Ji(7rg) - eJi(Treq)) 2 , (16) 

where the accessory parameter T = (1 — e 2 ) -2 . While vir- 
tually no difference within main spectral peak from the 
aperture filter of a circular aperture, the filter has a sig- 
nificantly greater transmission within the second and subse- 
quent peaks (e.g., for e = 0.3 about 2.5 times). As a result, 
the scintillation spectrum at high frequencies rises and the 
total scintillation power increases as well. 

Substituting the expression (|16p in the @, we obtain 
the weighting function Qs(z,e) which can be represented in 
the form: Qs{z, e) = Qs{z)Cs(e), where the integral Cs(e) is 
expressed in terms of hypergeometric functions as follows: 



Cs(e) = T 1 



15 v^r(|)r(f)£\(l,|;2;e 
3tt 



2 . 5/3 

■e +e ' 



(17) 
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Figure 5. Relative change in the scintillation power s 2 for SE 
regime (thick solid line, left scale) and for LE regime (dashed 
line, right scale) on the parameter e. The corresponding thin lines 
depict the approximations by the formulas 1 181 and 1201 
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Figure 6. Functions C(e, v) for e = 0,0.1,0.2,0.3,0.4,0.5,0.6 de- 
pending on v for SE (on left) and LE (on right) regimes. The solid 
lines depict VK model, dashed lines the GT model. 



There is a simple and good approximation for this function. 
Four members are enough to ensure an accuracy better than 
1% for e < 0.6: 

C s (e) « T(l + e 5/3 - 1.5682 e 2 - 0.1525 e 4 ) (18) 

The function Cs(e) is shown in Fig. [5] As can be seen, in the 
SE regime, the CO effect leads to a twofold increase in the 
power on the telescope with the e w 0.6. With widespread 
CO of 0.3 — 0.4, the effect is much smaller, though still of 
20 - 40%. 

In the LE regine, the CO effect is calculated similarly 
using the formula ([5]). The dependence of the scintillation 
on the parameter e is given by the formula 

c U \ - T (l 8 ![^W|£MzilM^ 2 , W\ 

c ^- T i i -27 — ifii — e+e )' 

(19) 

The simple approximation provides accuracy better than 
0.005 in the range < e < 0.6 can be written as 

C L (e) « T(l - 2.355 e 2 + e 8/3 + 0.2617 e 4 ) (20) 

The curve in Fig. [5] shows that in this regime, the CO does 
not introduce any significant effect for the most of telescopes. 
Even with the e = 0.6 the power increases only by 8%. The 
smallness of the effect is explained by the fact that for the 
Kolmogorov spectrum, there is some balance between the 
power decrease in the main peak and an increase in the 
second peak of the spectrum. 

Since the weighting function Qs(z,e) = Qs(z)Cs(e), all 
asymptotic dependencies on the altitude and diameter are 
retained. Accordingly, the effect on the measured scintilla- 
tion index s 2 is also described by the function C(e). It is 
worth noting that if the inner diameter of the aperture eD 
becomes comparable with the rp (telescope diameter is less 
than 0.5 m), then the presented above functions are approx- 
imate. 



4.2 The combined effect of the outer scale and 
the central obscuration 

The formulae, obtained in the previous section, are valid 
if the outer scale effect is negligible, i.e. dimensionless pa- 
rameter v < 1. Otherwise, the fraction of high-frequency 
scintillation spectrum grows and the CO effect increases as 
well. 

Obviously, after replacing the aperture filter in the ex- 
pressions for Qs(z,v) or Ql(z,v) (|12|l on (I16|l . the integral 
G{v) becomes dependent on the parameter e. Separation of it 
into two multiplicands, each depending on its own parame- 
ters, can not be performed. So, the combined impact of these 
effects can be written either as Q(z, e, v) = Q(z) Q(y, e)C(e) 
or as Q(z, e, v) = Q(z) Q(u)C(e, v). 

We chose the latter form, because for extremely large 
telescopes the effect of the outer scale is dominant. Accord- 
ingly the function C(e, v) depends on the adopted model and 
not only on the parameter v. The normalization, used in 
these calculations, provides that C(e, v) — C(e) when v = 0. 

The curves C(e, v), calculated for the VK and GT mod- 
els, are shown in Fig. [6] For both models the CO effect 
increases monotonically with v. It is seen, that for the VK 
model, which has a saturation at low spatial frequencies, 
the effect is more significant and at v ~ 1 almost reaches 
its maximum. This indicates that the scintillation within 
the main peak of the aperture filter, is suppressed almost 
completely. 

In the LE regime, the behavior of C(e,v) is also very 
interesting. For small v the CO effect is negligible on real 
telescopes, but at v ~ 0.5 becomes comparable with the 
effect in the SE regime and is growing further. Neverthe- 
less this does not compensate for the scintillation reduction 
caused by the outer scale effect (Fig. [4}. The VK and GT 
models differ more than in the SE regime. 
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Figure 7. Effect of the outer scale as a function of the diameter 
of telescope for two values of Lo for SE regime. The curve for 
the VK model is shown with the solid line, for GT model with 
the dashed, for EM with long dashed and for ME model with 
dash-dotted line. 



5 DISCUSSION 

5.1 Scintillation noise on large and and extremely 
large telescopes 

The value of the outer scale is important for many applica- 
tions, so it is often measured by different methods and its 
estimates and altitude profiles are available for s everal astro- 
nomi c al observatorie s |Abahamid et al. 2004 ; Maire et al. 
20071 ; iDali Ali et ail |20ld : iFlovd. Thomas-Osip fc Prietc 
20101 ). These data show that at altitudes greater than 8 km, 
the typical Lo is 20 — 25 m, although values from 10 to 100 m 
are sometimes observed. We will use these figures for further 
numerical evaluation. 

Fig.[7]and Fig.[8]show Q(y) as a function of the telescope 
diameter for Lo = 25 m and Lo = 50 m. Here the CO effect is 
not considered, so the curves are applicable only to telescope 
with e < 0.2. We see that even in SE regime on telescopes of 
10 m class, the scintillation power is « 0.6 from the power 
predicted by the standard relation (JSJ, assuming the actual 
spectrum of perturbations close to the VK or GT models. 
Exponential models lead to a smaller effect. 

For the designed extremely large telescopes TMT and 
E-ELT, the Q decreases to w 0.25, that certainly is to be 
taken into account in the estimates of errors budget in high- 
precision photometric observations. Note that the CO for 
the TMT (e = 0.11) and for the E-ELT (e = 0.15) is small, 
so that its accounting increases corrective functions by no 
more than 10 %. 

For the LE regime, the effect is more significant: the 
measured scintillations power is reduced by 10 — 20 times. 
Using the form ula ff]) and typical estimates of atmospheric 
moments from l|Kenvon et al.ll2006h . we can predict that for 
60 s exposure the scintillations noise becomes m 10 /imag 
instead of m 40 /imag predicted by the classical formula. For 
comparison, the similar estimate for 4 m telescope accounts 
~ 200 /imag. 

Naturally, the scintillation noise at the level of 10 /imag 
requires an appropriate photon noise. It is easy to calculate 
that 10 10 photons are to be accumulated at least, or, in the 
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Figure 8. Effect of the outer scale as a function of the diameter 
of telescope for two values of Lo for LE regime. See designations 
in Fig. [7] 



case of TMT, the source should provide illumination about 
25 photons cm - s _1 , which corresponds to a star with V « 
11.5 magnitude. Of course, for such photon fluxes, special 
technique should be used to avoid strong nonlinearity effects. 

Evaluation of the scintillation impact in the case of fast 
photometry (on time scales 0.001 — 0.01 s) for the same 
telescope, leads to the noise as 100 — 150 /imag, which is 
about 20 times less than for 4 m telescopes. The brightness 
of objects available for such measurements, depends on the 
exposure and for r = 0.01 s is equal to « 7.0 magnitude. 

The combined effect of the outer scale and the CO for 
two large (VLT and LSST) and two extremely large (TMT 
and E-ELT) telescopes is illustrated in Fig. [9l where the 
effect is represented as a function of the outer scale Lo. Ev- 
idently, the situation for the telescope LSST (D — 8.36 m) 
with a large e = 0.61 is distinguished. In the SE regime, the 
scintillation is larger than estimate gives by the formula ©, 
since CO effect dominates over the outer scale effect. How- 
ever, in the LE regime and for the most probable value of 
Lo — 25 m the scintillation power is expected half as many. 

The complex rarefied configuration of aperture of the 
GMT telescope has no evident diameter, so we can not ob- 
tain an estimate of the scintillation power in the usual way. 
Instead, we use the fact that the scintillation is almost uncor- 
related at dista nces greater than the diameter of a telescope 
(|Kornilovll2012ll . In the LE regime, the correlation can be as 
small as (~ 0.2), but owing to outer scale effect it should be 
further reduced. 

Ignoring the CO, although it for the central GMT mir- 
ror is large enough, we find that the scintillation power on 
GMT should be by 4 — 7 times smaller than for a single 
mirror with diameter 8.36 m. Therefore the curves for the 
telescope VLT, shown in Fig. O can be used as an estima- 
tion. 



5.2 General questions 

All the extremely large telescopes are designed with seg- 
mented primary mirror. The gaps between the individual 
segments and secondary mirror spiders have a characteris- 
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Figure 9. The combined effect of the outer scale and the CO as a 
function of L for 3 designed telescopes: LSST, TMT, E-ELT, and 
for the telescope VLT. The curves for the VK model are shown 
with solid lines, for GT model with dashed lines. The case of SE 
regime is shown on the left, the LE regime on the right. 



tic size of about 0.2 — 2 cm, and also increase the power 
of high-frequency scintillation. Maximum effect falls within 
second - third transmission peak of the aperture filter as 
in the case of the CO. In the first approximation, the effect 
is proportional to obscured relative area. The area of gaps 
between segments amounts to less than 1% (for TMT it is 
0.3%), so that even in the SE regime, we does not expect 
that their influence is noticeable. 

In assessing the observed effect of the outer scale in the 
section 13.21 we have assumed that the Lo is constant on 
the line of sight. Otherwise, the function Q(v} can not be 
factored out of integrals © or and this effect in the 
scintillation index should be evaluated for each layer and 
then integrate over the whole atmosphere. Evidently, it is 
necessary to know the vertical distribution of OT and Lo 
for this. 

However, the calculated effect is model-dependent and 
the imprecision, introduced by the uncertainty of the model, 
is quite large (see Fig. [3] and Fig. so variation of Lo with 
altitude on the order of ±50% can be neglected. The second 
circumstance, which has already been mentioned: in large, 
the effect is defined by the Lo in a bounded region of the 
altitudes (10 — 15 km) which generates the main part of the 
scintillations for large telescopes. 

The curves in Fig.[9]also explain that the temporal vari- 
ations of the measured scintillation power are to occur not 
only because of changes in OT intensity in the upper at- 
mosphere, but also because of variations in the outer scale. 
The net effect should be greatly depended on how these two 
parameters are connected. This question is virtually unex- 
plored, but it is clear that if an increase in the intensity 
correspond to a decrease of the Lo, then the scintillation 
power is to be quite stable characteristic, and vice versa. 

The functions, shown in Fig. [7] and [8] can be used for 
verification of type of the turbulence models. For this, the 
initial part of the dependencies (diameters of 2 — 8 m) for the 
SE regime is best suited. Here, simultaneous measurements 
with apertures of different diameters must be required. One 



can apply the known methods of integration of exit pupil 
image, selecting the desired sub-apertures during process- 
ing. In order to minimize an impact of the CO, this sub- 
apertures should be similar. E.g., on the VLT telescope one 
can distinguish a number of sub-apertures with a diameter 
of 8.2 m, 3.5 m and less with the same e = 0.15. 

Reduction of low-frequency components of the scintil- 
lation for non-Kolmogorov OT leads to yet another conse- 
qu ence: an effec tiveness of the pupil apodization, proposed 
m jYoundl 19671 ). should greatly increase. In the case of ex- 
tremely large telescopes, it is expected that optimal apodiza- 
tion should lead to a gain more than 2 times in the scintil- 
lation noise. 



6 CONCLUSION 

In this paper we considered two factors that modify the 
scintillation intensity in observations on large (4 — 10 m) 
and extremely large (20 — 50 m) telescopes. The first effect 
is caused by deviation of the real OT from the Kolmogorov 
model at the scales of the order of 10 m. In the models which 
describe the real turbulence with using outer scale Lo, the 
scintillation power at low spatial frequencies is much lower 
than for pure Kolmogorov spectrum. 

Effect of the outer scale can be described by an addi- 
tional function, which depends on the ratio of D/Lq. This 
function is equal to 1 when the telescope aperture is small 
and decreases with diameter increasing. We considered its 
general behavior by the example of piecewise power-law 
spectrum with the salient point at spatial frequency Lq -1 . 
Particular features were investigated by numerical integra- 
tion for four different models: Von Karman and Greenwood- 
Tarazano models, and two exponential models. 

For the observed values of Lo and models with a wide 
intermediate zone, the effect becomes visible on 4 m class 
telescopes. For measurements with long exposures this ef- 
fect is more important than in the case of short exposures, 
and for the extremely large telescopes TMT and E-ELT it 
can reduce the scintillation power ~ 10 times compared to 
classical estimates. 

For the Kolmogorov OT, the effect of amplification of 
the scintillation power, caused by the CO, inherent in every 
large telescope, results in the multiplication of the power 
for circular aperture by the function that depends only on 
the CO parameter. In the case of models with outer scale, 
the effect is described by a more complicated way, and it 
becomes significant for large apertures for both short and 
long exposures. 

The significant reduction of scintillation noise, due to 
the outer scale, enhances the potential of ground-based tele- 
scopes to study the variability of astronomical objects at 
~ 10 -5 level. This accuracy is sufficient to see, e.g., transit 
of Earth-like planet across the disk of Solar-like star. Drastic 
increase in the accuracy of fast photometry (up to ~ 10 -4 ) 
makes it possible to study the micro-variability of many as- 
tronomical objects without the accumulation of long time 
series suitable for temporal spectral analysis. 
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APPENDIX A: INTEGRALS 

Al Integrals in the case of Kolmogorov spectrum 

Here the known integrals, needed for calculating the weight- 
ing function, are listed with the formula ([4]) in SE regime: 



Is = 



POO 

/ q~ 11/3 q 3 (Mnq)) 2 d q 
Jo 



i 5 r(|)r3(|) 

and with formula ([5]) in LE regime: 



0.4508, (Al) 



r oo 

/ q- 11/3 q 2 (J 1 (nq)fdq 
Jo 

_ 27 2 2 / 3 V3r(f)r 2 (§) 
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= 0.8699. (A2) 



A2 Integrals in the case of finite limits 

In this appendix, the results of the integration over two spec- 
tral regions delimited by the dimensionless frequency v are 
written for the cases with a different exponent of the power 
spectrum of the phase perturbations <!>(?) fsee !3.2fl . Expres- 
sions for the SE and LE regimes are marked with the cor- 
responding subscript. The high-frequency part q > v with 
spectrum $(g) = q~ lx ^ 3 looks as 

r-oa 

Z3»= / q- 11/3 q 3 (Mnq)) 2 dq 

STY 2 7/3 , 7 3 13 „2,,2\ 

S_ ~28~ 3 ™' 2 ' T' 3 '-^ u ) ' 



I < , . = / q~ 11/3 q 2 {Ji(*q))'*l =" 8ir(l)r3(l: 



32vr 11 /6 

-^ 4/3 2 Fa(|,|i 1,2,3;-^). (A3) 

Low-frequency part q < v in the case of the spectrum 
with the saturation 4>(q) = v^ 11 ^ 3 : 



»= f ' v- 11 ' 3 q 3 (Ji(m))'- 

Jo 



dq 



,1/3 



((j^f + OA^)) 2 ), 



Xr>)= /V 1V V(Ji(7r<z)) 2 dg 
Jo 



= ^^ 3 2 f 3 (|,f;2,3,|;- 7 rV). (A4) 

Low-frequency part q < v in the case of the spectrum with- 
out the saturation $((?) = ^~ 2 g~ 5//3 : 

2s»= [" V- 2 q- 5/3 q 3 (Jl(Kq)) 2 dq 

Jo 
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= ^T- 4/3 ^3 (f,I; 2,3,1;-^). (A5) 



